Transverse Structure of Nucleon Parton Distributions from Lattice QCD 
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INTRODUCTION 

High energy lepton scattering reveals the quark and 
gluon structure of the nucleon by measuring matrix ele- 
ments of the light-cone operator 

O q (x) = J ^e iXx q{~n) ^Ve^^-^^^q^n). 

(1) 

The familiar quark distribution q(x) specifying the 
probability of finding a quark carrying a fraction x of the 
nucleon's momentum in the light-cone frame is measured 
by the diagonal nucleon matrix element, (P\0(x)\P) = 
q(x). Expanding Q(x) in local operators via the opera- 
tor product expansion generates the tower of twist-two 
operators, 

0{MiM2...M«} = ^>ii ^2 J ^M„l 5 t ( 2 ) 

and the diagonal matrix element l^p^Q^ 1 ^ 2 ---^^ spec- 
ifies the (n — l) th moment of the quark distribution 
/ dxx n ~~ x q{x). 

The generalized parton distributions H(x,£,t) and 
E(x, £, t) 0, Q, are measured by off-diagonal matrix 
elements of the light-cone operator 

(P'\0(x)\P) = ((A))H(x, £, t) + l -^((<7^n tl ))E(x, C, t), 

2m 

(3) 

where A^ = P'^ - P*, t = A 2 , | = -n ■ A/2, n 
is a light-cone vector, and ((r)) = U(P')TU(P). Off- 
diagonal matrix elements of the tower of twist-two opera- 
tors (P'|£)|' tl ' i, '"'' n ^|P) yield moments of the generalized 
parton distributions, which for £ = 0, are 

J dxx n - 1 H{x,Q,t) = A nQ {t) 

J dxx n ^E(x,0,t) = B n0 (t), (4) 



where A n i(t) and B n i(t) are referred to as generalized 
form factors (GFFs). 

Analogous expressions in which the light-cone operator 
O q {x) and twist- two operators contain an additional 75 
measure the longitudinal spin density, Aq{x) and spin- 
dependent generalized parton distributions H (x, £, t) and 
E(x, £, t) with moments A n i(t) and B n i(t). In this work, 
we present calculations of the generalized form factors 
A/ n —i t 2,s),o(t) and ^4(n=i.2,3),o(*) m full QCD and discuss 
their physical significance. 



TRANSVERSE STRUCTURE OF PARTON 
DISTRIBUTION 

In general, H (x, £, t) is complicated to interpret phys- 
ically because it combines features of both parton distri- 
butions and form factors, and depends on three kinemat- 
ical variables: the momentum fraction, x, the longitu- 
dinal component of the momentum transfer, £, and the 
total momentum transfer squared, t. In the particular 
case in which £ = 0, however, Burkardt Q has shown 
that H(x,0,t), as well as its spin-dependent counterpart 
H(x,0,t), has a simple and revealing physical interpre- 
tation. 

It is useful to consider a mixed representation in which 
transverse coordinates are specified in coordinate space, 
the longitudinal coordinate is specified in momentum 
space, and one uses light-cone coordinates for the lon- 
gitudinal and time directions: x^ — (x° ± x 3 )/^/2, 
p± _ ±p 3 y,/2. Using these variables, letting x de- 
note the momentum fraction and b± denote the trans- 
verse displacement (or impact parameter) of the light- 
cone operator relative to the proton state, one may de- 
fine an impact parameter dependent parton distribution 
in light-cone gauge 

q(x,b x ) = (P + ,Rx=0,X\O q (x,b x )\P+,R x =0,X), (5) 
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where 

o q (xj ± )= J ^ xp+x ~q(-^Ah + q (%r,b ± )- (6) 

Burkardt shows that the generalized parton distribu- 
tion H(x,0,t) is the Fourier transform of the impact pa- 
rameter dependent parton distribution, so that 



ff(a;,0,-Ai) = / d 2 b ± q(x,b ± )e lb±A± , 

A n0 (-Al) = f d 2 b ± [ dxx^qixj^e^^, (7) 



where the second form follows from Eq. ©. Although 
one normally only expects a form factor to reduce to 
a Fourier transform of a density in the non-relativistic 
limit, Ref. 4] shows that special features of the light- 
cone frame also produce this simple result in relativistic 
field theory. Thus, H (x, 0, t) specifies how the transverse 
distribution of quarks varies with the longitudinal mo- 
mentum fraction x. 

Physically, we expect the transverse size of the nucleon 
to depend significantly on x. Averaging q(x, b±) over all 
x produces Aifi(t) and thus corresponds to calculating 
the form factor. Hence, the average size is character- 
ized by the transverse rms radius (rj_) 5 = (x\ + x 2 }? = 

■^/|(r 2 )ii. From the experimental electromagnetic form 
factor, the transverse rms charge radius of the proton is 
0.72 fm. As x — * 1, a single active parton carries all the 
momentum and the spectator partons give a negligible 
contribution. In this case the active parton represents the 
(transverse) center of momentum, and the distribution 
in impact parameter reduces to a delta function S 2 (b±) 
with zero spatial extent. Indeed, explicit light-cone wave 
functions jjj Q bear out this expectation, with the result 



7] 



q(x,bj 



n.c a 



Y\_ dxj d 2 

3 = 1 



±3 




fj_j \ 6 (x - x a ) 



x)r± a 



X *n,c(^i: ■ ■ -;^xij • ■ -)*n,c(^i, ■ • • . . .), 

where a denotes the index of the active parton, n is the 
number of partons in the Fock state and the sum over 
c represents the sum over all additional quantum num- 
bers characterizing the Fock state. Here, one explicitly 
observes lim^i q(x, b±) cc S 2 (b±). Since H(x, 0, t) is the 
Fourier transform of the transverse distribution, the slope 
in —t = at the origin measures the rms transverse 



radius. As a result, we expect the substantial change in 
transverse size with x to be reflected in an equally sig- 
nificant change in slope with x. In particular, as x — > 1 
the slope should approach zero. Hence, when we calcu- 
late moments of H (x, 0, £), the higher the power of x, the 
more strongly large x is weighted, and the smaller the 
slope should become. Therefore, this argument makes 
the qualitative prediction that the slope of the general- 
ized form factors A n0 (t) and A n0 (t) should decrease with 
increasing n, and we expect that this effect should be 
strong enough to be clearly visible in lattice calculations 
of these form factors. 



LATTICE MEASUREMENT 

We consider the three spin independent moments 
(P'\0^\P) = (( 7 w Hio(t) 



(P'\0 { ^ 2} \P) 



— <(<7^)>A Q B 10 (t) , 



P { ^((r 2} ))A 20 (t) 

2-p{^((a^ a ))A a B 2 o(t) 
2m 



+ 2-p{^P^((a^ a ))A a B 30 (t) 
2m 

+ A^A^(( 7 ^>))A 32 (i) 

+ ^-A^W({a^ a ))A a B 32 (t), (8) 
2m 

where P M = (P^ + P u )/2, as well as the analogous spin 
dependent moments |8j. 

Generalized form factors ^(n=i,2.3),o(0 an d 
A(n=i,2,3),o(t) were calculated using the new method in- 
troduced in Ref. |9( . We considered all the combinations 
of P and P' that produce the same four-momentum 
transfer t = (P' — P) 2 , subject to the conditions that 
P = mS n ^n y ,n z ) and P' = (0,0,0) or ^-(-1,0,0). 
Using all these momentum combinations for a given t 
below 3.5 GeV, we calculated all the hypercubic lattice 
operators and index combinations that produce the 
same continuum GFFs, obtaining an overdetermined 
set of equations from which we extracted a statistically 
accurate measurement. The errors are substantially 
smaller than obtained by the common practice of 
measuring a single operator with a single momentum 
combination. Our calculations are based on unqucnchcd 
SESAM configurations 10] on 16 3 x 32 lattices with 
k = 0.1560 and k = 0.1570, corresponding to pion 
masses of m^ = 897 and 744 MeV respectively. 
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FIG. 1: Normalized lattice results for generalized form factors 
A„o and A„q as a function of momentum transfer squared, —t, 
for n=l (circles), 2 (triangles), and 3 (squares). 



and A 2j o depend on the sum. Hence only moments dif- 
fering by two compare the same physical quantity with 
different weighting in x. To facilitate determination of 
the slope of the form factors and to guide the eye, the 
data have been fit using a dipole form factor 



.4 



dipole 



.4 



1 - 



(9) 



The solid line denotes the least-squares fit and the shaded 
error band shows the error arising from the statistical 
error in the fit mass, Am^. Although the dipole fit is 
purely phenomenological, we note that it is consistent 
with the lattice data. For reference, the normalization 
factors A u q and dipole masses are tabulated in Table HJ 

The top panel in Fig. 1 shows the flavor non-singlet 
case A u — A d , for which the connected diagrams we have 
calculated yield the complete answer. It is calculated at 
the heaviest quark mass we have considered, correspond- 
ing to = 897 MeV. Note that the form factors are 
statistically very well separated, and differ dramatically 
for the three moments. Indeed, the slope at the origin 
decreases by more than a factor of 2 between n = 1 and 
ii = 3, indicating that the transverse size decreases by 
more than a factor of 2. The second panel shows analo- 
gous results for lighter quarks, — 744 MeV, where we 
observe the same qualitative behavior but slightly weaker 
dependence on the moment. The third panel shows the 
flavor singlet combination A u +A d , for which we have had 
to omit the disconnected diagram because of its signifi- 
cantly greater computational cost. Comparing this figure 
with the top panel calculated at the same quark mass, we 
observe that while the connected contributions to A u ±A d 
are qualitatively similar, there is significant quark flavor 
dependence that can be used to explore the nucleon wave 
function. It is useful to note our results for the u and d 
GFFs are consistent with the n=2 moments calculated 
in Ref. The bottom panel shows the spin-dependent 
flavor non-singlet form factors A u — A d at the heaviest 
quark mass. Thus, comparing the top and bottom fig- 
ures displays the difference between the spin averaged 
and spin dependent densities. We observe a striking dif- 
ference, in that the change between the n = 1 and n = 3 
form factors for q(x,b±)^ — q(x,b±)i is roughly 6 times 
smaller than for ^(q(x, + q(x, b±)i). 

Finally, it is useful to use the slope of the form factors 
at t = to determine the transverse rms radius, 



Figure 1 presents our principal results, showing the 
generalized form factors A n o(t) and A n o(t) for the low- 
est three moments: n = 1, 2, and 3. The form factors 
have been normalized to unity at t = to make the de- 
pendence of the shape on n more obvious. Note that 
^4-1,0:^3,0: an d ^.2,0 depend on the difference between 
the quark and antiquark distributions whereas ^4i,o, ^4.3 0) 



2 



(ri) 



J^b^blJdxx^qix,^) 
Jd 2 bj_ J dxx n - 1 q(x,b±) 



(10) 



Transverse rms radii calculated in this way are tabulated 
in Table [i] To set the scale, the transverse charge ra- 
dius at this mass is (fj^) charge — 0-48 fm, which is two- 
thirds the experimental transverse size 0.72 fm, reflecting 
the absence of a significant pion cloud. For the heaviest 
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GFF 


A(0) 


m d (GeV) 


(r'i)i ( fm ) 


= 897 MeV (k = 0.1560) 


A i ;- d (o) 

1,1) V / 


1.000 ± .001 


1.470 ± .031 


0.380 ± .008 


A?;- n d (0) 


0.241 ± .004 


2.102 ± .081 


0.266 ± .010 


A^ d (Q) 


0.060 ± .008 


3.857 ± .494 


0.145 ± .019 




2.998 ± .002 


1.205 ± .014 


0.463 ± .005 


A-+ a (0) 


0.666 ± .009 


1.706 ± .040 


0.327 ± .008 




0.155 ± .018 


2.099 ± .153 


0.266 ± .019 


A^ o d (0) 


1.195 ± .014 


1.850 ± .028 


0.302 ± .005 


A-- O d (0) 


0.293 ± .006 


2.223 ± .058 


0.251 ± .007 


^3,o d (0) 


0.123 ± .004 


2.233 ± .087 


0.250 ± .010 


= 744 MeV (k = 0.1570) 


^,o"(0) 


1.001 ± .001 


1.402 ± .019 


0.398 ± .005 


A^ d (0) 


0.261 ± .009 


1.814 ± .049 


0.308 ± .008 


^3,o"(0) 


0.071 ± .013 


2.373 ± .138 


0.235 ± .014 



TABLE I: Generalized form factors at t = 0, dipole masses, 
and transverse rms radii for the cases plotted in Fig. 1. 

quark mass, m n — 897 MeV, the non-singlet transverse 
size (^) u -(j = 0.38 fm is slightly smaller than the rms 
charge radius, but drops 62% to 0.14 fm for n—3. The 
singlet size (r^) u+c j is 0.46 fm, and drops 43% to 0.27 for 
ri=3. This is a truly dramatic change in rms radius aris- 
ing from changing the weighting by x 2 . An alternative 
way to describe the same effect is in terms of the mean 
value of x. The mean value of x in the distribution q(x) 
is of the order of 0.2 and roughly 0.4 in the distribution 
x 2 q(x). In these terms, the non-singlet transverse size 
drops 62% as the mean value of x increases from 0.2 to 
0.4, and goes to zero when x reaches 1. 

CONCLUSIONS AND OUTLOOK 

In the "heavy pion world" presently accessible to full 
lattice QCD, we have calculated the lowest 3 generalized 
form factors A n o and A n o up to \t\ = 3 GeV as shown 
in Fig. H We obtain excellent precision for n = 1 and 
sufficient precision for n — 2 and 3 to clearly distinguish 
each form factor and observe striking differences in slope 
and hence transverse size. Whereas there are other cal- 
culations of isolated moments, three moments are crucial 
for the present investigation since n= 1 and 3 are neces- 
sary to measure the same combination of quark and an- 
tiquark distributions. The dependence of the transverse 
size on x is most dramatic for the heaviest u — d combi- 
nation, for which {r\) n _^ decreases by 62% between the 
first and third moment. We also observed clear depen- 



dence of the transverse distribution on flavor and spin. 
Our results show that the commonly used factorization 
Ansatz H(x,0,t) = Q(x)F(t) is fundamentally wrong in 
the "heavy pion world" and we are aware of no arguments 
as to why it should be restored for lighter quarks. 

The most immediate challenges are to extend these cal- 
culations to the chiral regime of realistic quark masses, 
which is being explored using a hybrid calculation of dy- 
namical staggered sea quarks and domain wall valence 
quarks [12|, and to extend techniques for evaluating dis- 
connected diagrams 0] to GFFs. When precise, con- 
trolled extrapolations to the physical pion mass are fi- 
nally achieved, moments calculated from first principles 
will play an essential role in complementing experimental 
results because of the impracticality of measuring the full 
x, £, and t dependence of H(x,^, t) and H(x,£,,t) exper- 
imentally. In addition, they will provide rich insight into 
the flavor and spin dependence of the transverse wave 
function. 
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